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Spectroscopy of higher bottomonia 
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In this contribution, we discuss our recent unquenched quark model results for the spectrum of 
bb mesons with self energy corrections, due to the coupling to the meson-meson continuum. Our 
unquenched quark model predictions for the masses of the recently discovered Xb{3P) states are 
compared to those of a re-fit of Godfrey and Isgur’s relativized quark model to the most recent 
experimental data. The possible importance of continuum effects in the \b(3P) states is discussed. 
Finally, we show our quark model results for the radiative decays of the Xb(3P) system and the 
open-bottom decays of bb mesons. 


I. INTRODUCTION. AN UNQUENCHED 
QUARK MODEL FOR MESONS 


between these low-lying states and the first open-bottom 
decay thresholds. 


In the UQM for hadrons [MI , the effects of qq sea pairs 
are introduced explicitly into the quark model (QM) via 
a QCD-inspired 3 Po pair-creation mechanism |8(. The 
qq pairs have vacuum quantum numbers and their effect 
on hadron observables can be added as a perturbation. 
This approach is a generalization of the unitarized quark 
model by Tornqvist and Zenczykowski 0 and was also 
motivated by the work of Geiger and Isgur la¬ 
in the UQM for mesons 0-0 , the meson wave function, 
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is written as the superposition of a zeroth order quark- 
antiquark configuration and a sum over the possible 
higher Fock components, due to the creation of 3 Po qq 
pairs. Above threshold, this coupling to the continuum 
leads to open-flavor strong decays; below threshold, it 
leads to virtual qq — qq components in the meson wave 
function and a shift of the physical meson mass with re¬ 
spect to its bare mass, namely a self-energy correction. 
Something similar was also done in the baryon sector jlj- 
0], with the intention to extend the ran ge of ap plicability 
of the baryon QM in its various forms [Ill-l20l|. 

In Ref. 0 ], the UQM formalism was used to calculate 
the mass shifts of IS 1 , 2 S and IP bottomonium states due 
to the coupling to the meson-meson continuum. This was 
just a preliminary study, in which the bare energies were 
taken as free parameters and their sum with the self- 
energies fitted to the physical masses of the mesons of 
interest. There, we also calculated the probability to find 
the wave function of IS 1 , 2 S and IP bottomonium states 
in the continuum component. The results vary between 
2% and 6% and their smallness is due to the distance 
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In Ref. 0, the UQM was used to calculate the char- 
monium spectrum (IS 1 , 2 S, IP, 2P and ID states) with 
self-energy corrections. Unlike Ref. 0 ], the bare energies 
were calculated explicitly within a potential model 1211 ] . 
thus extending the predictive power of the UQM. Our re¬ 
sults were used to discuss the nature of the A (3872) [20 ] . 
a meson characterized by J pc = 1 ++ quantum numbers 
HU, a very narrow width (less than 1.2 MeV HU), and a 
mass 50 — 100 MeV lower than QM predictions. Accord¬ 
ing to our study, the X(3872) can be interpreted as a cc 
core plus higher Fock components due to the coupling to 
the meson-meson continuum. The problem of the mass 
difference between QM predictions HI, !H 0 and the 
experimental data for the A(3872) was solved thanks to 
a downward mass shift, due to self-energy corrections. 
See also Ref. 0. 


In Ref. 0, this calculation was extended to bottomo¬ 
nia and used to provide predictions for IS, 2 S, 3 S, IP, 
2P, 3P and ID resonances; the open-bottom decays of 
bb states were calculated too within the 3 Po pair-creation 
model 0. Analogously to what happens for the X(3872) 
in the cc sector, threshold effects may play an impor¬ 
tant role in the Xb(3P) system, discovered by ATLAS in 
2012 H3, HU: because Xb(3P) mesons are very close to 
the first open-bottom decay thresholds. Thus, to distin¬ 
guish between QM and UQM interpretations for these 
mesons, our UQM results for the mass barycenter of the 
Xb(3P) system were compared to those of a re-fit of God¬ 
frey and Isgur’s model mass formula to the most recent 
experimental data. In Ref. 0, the radiative transitions 
of Xb(3P) mesons were calculated too. 


In this contribution, we briefly analyze the results of 
Refs. 00 for the bottomonium spectrum, the open- 
bottom decays and the radiative transitions of Xb(3P) 
mesons. These results are then used to discuss the pos¬ 
sible importance of continuum effects in Xb(3P) mesons. 
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FIG. 1: Spectrum of bb mesons in the UQM. Picture from 
Ref. [?|; APS copyright. The calculated masses of bb states 
are shown as black lines, the experimental data [24J as green 
boxes. 


II. SPECTROSCOPY OF HIGHER 
BOTTOMONIA 


A. Bottomonium spectrum in the UQM 


Here, we discuss our UQM results for the spectrum of 
bb mesons with self energy corrections of Refs. 00 - In 
our calculation, we fitted 


M a =E a + S(E a ) (2) 

to the experimental data (24[. The physical mass of a 
meson, M a , results from the sum of a bare energy term, 
Eg, computed within Godfrey and Isgur’s relativized QM 
I2ll |. and a self energy correction, 


£(K) = E / ^ 


BC ' 
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evaluated in the unquenched quark model formalism 
0,11-0. See Fig. |T] In Ref. 0, we also gave results 
obtained by re-fitting the relativized QM mass formula 
[2l| | to the most recent experimental data jH], without 
including continuum effects. The re-fit was necessary to 
compute the open-bottom strong decays, which required 
more precise values for the masses of the higher-lying 
decaying mesons, like 45 states. 

It is interesting to notice that (see Refs. 00): ^be¬ 
low 10.5 GeV, our UQM and relativized QM results are 
very similar; 2) above 10.5 GeV, an evident discrepancy 
emerges in the case of x'fc(3P) resonances (20 [28|, due to 
the vicinity of these states to the first open-bottom decay 
thresholds. This means that, far from thresholds, contin¬ 
uum effects can be (at least partially) re-absorbed into 
a redefinition of the model parameters, as also demon¬ 
strated in Ref. [29|; close to thresholds, this is not possi¬ 
ble anymore and continuum effects can play an important 
role. For exa mple, t his also happens in the case of the 
V(3872) 00|3(Hi! and 19^(2317) [H mesons. 


TABLE I: Mass barycenters of Xf,(lP)> xt(2P) and x&(3 P) 
systems (column 1 ) and mass splittings between the mem¬ 
bers of the x&(lP), Xb{2P) and yt(3 P) multiplets (column 2 
and 3) in the UQM. See Ref. ]7j. The results are expressed in 
MeV. The notation AAf 2 i(lP) stands for the mass difference 
between the x& 2 (l P) and xm( 1R) resonances, AMio(lP) for 
the mass difference between the Xw(lU) and Xw(lP) reso¬ 
nances, and so on. 


M Xb(1P ) AM 2 i( 1P) AMro(lP) 
9894 21_30_ 

Af X6 ( 2 p) AM2i(2P) AA-/io(2P) 
10241 21_ 14 _ 

M Xb(3P) AA/ 2 i(lP) AAfio(lP) 
10510 17 13 


TABLE II: Mass barycenters of Xf>(lP), Xt>(2P) and Xfc(3P) 
systems (column 1 ) and mass splittings between the members 
of the X&(1U), Xb{2P) and xt(3P) multiplets (column 2 and 
3), from our re-fit of Godfrey and Isgur’s mass formula of Ref. 
[ 3 . The results are expressed in MeV. 


B. Mass barycenter of the Xf>(3P) system 

Our UQM results were also used to evaluate the mass 
barycenter of the Xb(3P) system. See Table HI An¬ 
other estimation was obtained by re-fitting the rela¬ 
tivized QM mass formula [2l| to the most recent ex¬ 
perimental data [H, [ 34 ]. See Table [TT] The two sets 
of theoretical results can be then compared to the ex¬ 
perimental data: M Xb ^ P ^ = 10.534 ± 0.009 GeV 24], 
M Xh(3P) = 10.530 ± 0.005(stat.) ±0.009(syst.) GeV [27f, 
and M Xb{3P) = 10.551 ± 0.014(stat.) ± 0.017(syst.) GeV 
[HI]. LHCb has also recently provided a measure for 
the mass of the Xbi(3P) resonance [35[: M Xbl ( 3 p\ = 
10515.7^3;g(stat.)l2'i( s yst-) MeV. See also Ref. 36l] , 
where the authors provided predictions for the splittings 
between the members of the xt,(3P) multiplet using sev¬ 
eral potential models. 


A7 Xc (i p) AM 21 UP) AMio(lP) 
3494 46_96_ 

Af Xc (2 P ) AA42i(2P) AA/io(2P) 
3906 56 -47 


TABLE III: Experimental data [HJ for the mass barycenters 
of Xc(lP) and Xc(2P) systems (column 1) and mass splittings 
between the members of the Xc(lR) an d Xc(2P) multiplets 
(column 2 and 3). The results are expressed in MeV. 
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Meson 

Mass [MeV] 

J^ 

BB 

BB* 

/ BB* B*B* 

B S B S 

B S B* 

/ B S B* S B* a B* s 

T( 10580) or 

T(4 il 5'i) 10.595 

10579.4 ± 1.2 f 

1— 

20 

— 

— 

— 

— 

— 

Xb2(2 3 P 2 ) 

10585 

2 ++ 

34 

- 

- 

- 

- 

- 

T(3 3 Di) 

10661 


23 

4 

15 

- 

- 

- 

T 2 (3 3 D 2 ) 

10667 

2 "” 

- 

37 

30 

- 

- 

- 

r 2 (3 1 D 2 ) 

10668 

2 "+ 

- 

55 

57 

- 

- 

- 

T 3 (3 3 D 3 ) 

10673 

3" 

15 

56 

113 

- 

- 

- 

Xw(4 3 Po) 

10726 

0 ++ 

26 

- 

24 

- 

- 

- 

T 3 (2 3 G 3 ) 

10727 

3" 

3 

43 

39 

- 

- 

- 

Xw(4 3 Pi) 

10740 

1 ++ 

- 

20 

1 

- 

- 

- 

/ii,(4 1 Pi) 

10744 

1 +- 

- 

33 

5 

- 

- 

- 

X62(4 3 P 2 ) 

10751 

2 ++ 

10 

28 

5 

1 

- 

- 

Xb2(3 3 P 2 ) 

10800 

2 ++ 

5 

26 

53 

2 

2 

- 

X 3 (3 1 Ph) 

10803 

3+^ 

- 

28 

46 

- 

3 

- 

T(10860) or 

T(5 3 Si) 10876 ±ll f 

1 " 

1 

21 

45 

0 

3 

1 

T 2 (4 3 D 2 ) 

10876 

2 

- 

28 

36 

- 

4 

4 

T 2 {4 1 D 2 ) 

10877 

2~+ 

- 

22 

37 

- 

4 

3 

T 3 (4 3 D 3 ) 

10881 

3- 

1 

4 

49 

0 

1 

2 

T 3 (3 3 G 3 ) 

10926 

3— 

7 

0 

13 

2 

0 

5 

T(11020) or 

T( 6 3 S'i) 11019 ± 8 f 

1~~ 

0 

8 

26 

0 

0 

2 


TABLE IV: Open-bottom decay widths (in MeV) for higher bottomonium states, as from Ref. [7}; APS copyright. Column 2 
shows the values of the masses of the decaying bb states: when available, we used the experimental values from PDG [24| (f), 
otherwise the theoretical predictions of the relativized QM [21,] , whose mass formula we re-fitted to the most recent experimental 
data 01 • Columns 3-8 show the decay width contributions from various BC channels, like BB, BB* and so on. The values of 
the 3 P 0 model parameters were fitted to experimental data for the strong decays of bb resonances 0. The symbol - in the table 
means that a certain decay is forbidden by selection rules or that the decay cannot take place because it is below the threshold. 


It is worthwhile noting that: 1) our UQM and rel¬ 
ativized QM predictions of Tables Q] and eh are simi¬ 
lar in the case of the mass barycenter, which is the 
average between the three meson masses of the mul- 
tiplet, but differ for the splittings between the mem¬ 
bers of the multiplet; 2) the mass splitting scheme, 
with M Xbl ( 1P2 p) — M Xb0 (ip 2 p) « 15 — 30 MeV and 
M Xb2 (iP 2 P) — M Xbl ( 1Pt2 p) ~ 15 — 20 MeV, in the UQM 
is broken for Xb(3P) mesons. Something similar also 
happens for x c (2P) mesons. See Table IIIII Therefore, 
we think that new and more precise experimental data 
will be necessary to confirm or discard our hypothesis of 
Xb(3 P) mesons as made up of a bb core plus non-negligible 
continuum components, due to threshold effects. 


III. OPEN BOTTOM STRONG DECAYS IN 
THE 3 P 0 PAIR-CREATION MODEL 

Here, we discuss our results for the open-bottom strong 
decays of Ref. [tJ . The transitions we considered are of 
the type 


A -a B + C , . 

bb bq qb ’ 1 J 

where q is a light quark (u, d or s). The decay widths 
were computed as [37[ : 

r a^bc = &A^>Bc(qo) E \{BCq 0 U\T^\A)\ 2 . (5a) 
i,j 


Here, T f is the pair-creation operator of the 3 Pq model, 
where we substituted the constant pair-creation strength, 

rr 

7 o, with an effective one, Tn", to suppress unphysical 
heavy quark pair-creation l38j |, and introduced a 
Gaussian quark form factor, because the pair of created 
quarks has an effective size Rl }?j. [Tol l29i. ••{!) . ^A-^Rn(qn) 
is the standard relativistic phase space factor 0-0, [37J , 


‘f’T-s-BC 


27t<7o 


Eb(qo)E c (q 0 ) 
M a 


(5b) 


which depends on the relative momentum go between B 
and C and on the energies of the two intermediate state 
mesons, E b = M b + q% and E c = \JMj? + q%. The 
values of the masses M a , Mb and M r , used in the cal¬ 
culation, were taken from the PDG [24j and Ref. ]34j . 
while, in the case of still unobserved states, we used the 
predictions from our re-fit of Godfrey and Isgur’s model 
mass formula of Ref. 7]. 


IV. QM CALCULATION OF x 6 (3P)’S El 
RADIATIVE TRANSITIONS 

In this section, we show our QM results for the radia¬ 
tive transitions of Xh(3P) states from Ref. [g|. The El 
radiative transitions were calculated according to [25] : 

4 

Tei = -C fz S S s'e 2 b a\^ f \r\^)\ 2 E 3 ^- . 

3 M} bb) 


(6a) 
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Transition _ P 7 (QM) (MeV) r b5 (QM) (KeV) P 7 (UQM) (MeV) T^(UQM) (KeV) 


X6o(3 j P 0 ) -4 T(1 j Si)7 983 

0.6 

984 

0.2 

X6o(3 3 P 0 ) -4 T(2 3 Si ) 7 460 

1.2 

461 

0.7 

X6o(3 3 Po) -4- T(3 3 Si ) 7 138 

6.1 

139 

6.4 

X6o(3 3 P 0 ) -4 T( 1 3 Di ) 7 344 

0.2 

376 

0.2 

X6o(3 3 Po) —> T( 2 3 Di ) 7 69 

0.9 

- 

- 

X6i(3 3 Pi)^T(l 3 Si ) 7 971 

2.1 

1061 

1.7 

X6i(3 3 Pi) T(2 3 Si)7 477 

2.5 

542 

2.7 

X6i(3 3 Pi) -4 T(3 3 Si ) 7 155 

7.4 

223 

22.8 

X&i(3 3 Pi) -4- T(1 3 P>i ) 7 361 

0 

458 

0 

X6i(3 3 Pi) -4- T( 2 3 Ph )7 86 

0.4 

- 

- 

Xbi(3 3 Pi) -4- T 2 ( 1 3 P 2 )7 341 

0 

408 

0.1 

Xm(3 3 Pi) ->• T 2 (2 3 D 2 ) 7 79 

1.0 

- 

- 

X6 2 (3 3 P 2 ) -4- T(1 3 Si)7 1010 

3.9 

1059 

3.6 

Xt2(3 3 P 2 ) —>• T(2 3 Si ) 7 489 

3.8 

540 

4.0 

X6 2 (3 3 P 2 ) -4 T(3 3 Si ) 7 168 

8.2 

221 

19.6 

X6 2 (3 3 P 2 ) T(1 3 Di ) 7 373 

0 

456 

0 

X6 2 (3 3 P 2 ) -4- T(2 3 P>i ) 7 99 

0 

- 

- 

X62(3 3 P 2 ) -4- T 2 (1 3 D 2 )7 354 

0 

406 

0 

X6 2 (3 3 P 2 ) -4- T 2 (2 3 P 2 )7 92 

0.3 

- 

- 

X62(3 3 P 2 ) -4- T 3 (1 3 P 3 )7 357 

0 

441 

0 

X6 2 (3 3 P 2 ) -4- T 3 (2 3 P 3 )7 86 

1.4 

- 

- 


TABLE V: PI radiative transitions of Xb(3P) states, calculated with Eq. © [t|. In the calculation, we used our QM and 
UQM predictions for the masses of the \b{ 3P)’s fj] and the wave functions of the relativized QM [2l| . Table from Ref. Q; 
APS copyright. 


Here, et, = — ^ is the charge of the 6-quark (in units of 
e), a is the fine structure constant, E 1 is the final photon 

energy, is the total energy of the final 66 state, M'j b ^ 
is the mass of the initial 66 state, (i/’/l r \ipi) is a radial 
matrix element, which involves the initial and final radial 
wave functions, and the angular matrix element Cfi is 
given by 

C fi = max(P, L')(2J' + 1) | ^ ^ f J . (6b) 

We calculated the matrix elements of Eqs. (© assum¬ 
ing, for the initial and final states, the wave functions of 
Godfrey and Isgur’s relativized QM HU], We gave results 
considering, for the masses of the decaying \{,(3P) reso¬ 
nances, predictions from the UQM and our re-fit of the 
relativized QM mass formula of Ref. [7). This was done 
because the phase space factor of the decay strongly de¬ 
pends on the masses of the initial and final state mesons 
and our UQM and relativized QM predictions for the 
masses of Xb(3P) states are quite different. 

Finally, our results are reported in Table l~Vl 

V. CONCLUSION 

In this contribution, we discussed the results of an 
UQM calculation of the bottomonium spectrum with self¬ 


energy corrections 0 . In the UQM 00 , the effects of qq 
sea pairs are introduced explicitly into the quark model 
(QM) via a QCD-inspired 3 Po pair-creation mechanism 
[S|. The qq pairs have vacuum quantum numbers and 
their effect on hadron observables can be added as a per¬ 
turbation. Our UQM results are similar to those of the 
relativized QM [2lJ for the majority of states we consid¬ 
ered, but differ substantially in the case Xb(3P) mesons. 
This happens because \b(3P) states are very close to the 
first open-bottom decay thresholds, thus meson-meson 
continuum effects can play an important role. A com¬ 
parison between UQM and relativized QM results for 
the mass barycenter of Xfc(3P) states was done but, up 
to now, the present experimental data do not permit to 
conclude which of the two interpretations is preferable. 
Finally, we discussed our results for the open-bottom 
strong decay widths of 66 states within the 3 Po pair- 
creation model 0 and for the radiative transitions of 
Xb(3P) states in the QM formalism Q. 


In Refs. H3, the QM was extended the include the 
quark-antiquark-gluon higher Fock components to study 
the spectroscopy of hybrid mesons. In Ref. HU can be 
found a classification of possible tetraquark states. 
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10.5 - 

10 . 0 - 
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0 -+ 1— 1+- o ++ 1++ 2++ 2— 







